Spin relaxation in CdTe quantum dots with a single Mn atom 
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We have investigated spin relaxation times in CdTe quantum dots doped with a single Mn atom, 
a prototype of a system where the interaction between a single charge carrier and a single spin 
takes place. A theoretical model that was used includes the electron-Mn spin exchange interaction 
responsible for mixing of the states of different spin in the basic Hamiltonian and electron-phonon 
interaction as a perturbation responsible for transitions between the states. It was found that the 
dominant electron-phonon interaction mechanism responsible for spin relaxation is the interaction 
with acoustic phonons through deformation potential. Electron and Mn spin relaxation times at 
room temperature take values in the range from microseconds at a magnetic field of 0.5 T down 
to nanoseconds at a magnetic field of 10 T and become three orders of magnitude larger at cryo- 
genic temperatures. It was found that electron spin-orbit interaction has a negligible effect on spin 
relaxation times, while the changes in the position of the Mn atom within the dot and in the dot 
dimensions can change the spin relaxation times by up to one order of magnitude. 
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I. INTRODUCTION 

The potential of utilizing the spin degree of freedom as 
a classical bit in spintronic devices^— and as a quantum 
bit in potential quantum information processing devices 4 - 
has been largely recognized in the last two decades. Op- 
eration of these devices crucially depends on the ability 
to manipulate the spin degree of freedom of the system^ 
However, in all realistic open systems, undesirable spin 
flips occur due to the interaction with environment. It 
is therefore of great importance to understand and be 
able to quantitatively describe the mechanisms of spin 
dephasing and spin relaxation. 

Physical systems which are expected to be particularly 
suitable for aforementioned applications are based on ar- 
chitectures that contain quantum dots - artificial nanos- 
tructures where charge carriers are confined in all three 
spatial directions i&l Due to quantum confinement effect 
the spectrum of electronic states in quantum dots is dis- 
crete and as a consequence phase space for relaxation and 
dephasing processes is greatly reduced^ Therefore, long 
spin lifetimes of electrons confined in quantum dots are 
expected. 

In addition, quantum dots provide a playground where 
one can study fundamental interactions on a single carrier 
or spin level. Such a level of understanding is necessary 
before one can proceed to understand more complicated 
device structures. Quantum dots doped with a single 
Mn atom have drawn particular attention for fundamen- 
tal studies in recent years. The manganese atom acts 
effectively as an additional spin 5/2 degree of freedom 
and therefore enables fundamental studies of interaction 
between a single charge and a single spin in these dots. 
Experimentally observed signature of this interaction is 
the splitting of an cxciton line in quantum dot photolu- 
minescence spectrum^— Predictions that it is possible 
to optically manipulate the state of a Mn spin in the 
quantum do t 12 ' 14 ' 15 were recently realized in several ex- 
periments*^— 

However, very little is known about the lifetimes of Mn 



and electron spin in these dots. A theoretical analysis 1 ^ 
of optical orientation experiments 1 ^ provided estimates 
of hole and electron spin relaxation times necessary for 
Mn spin orientation to occur. In a very recent work, Mn 
spin relaxation times in quantum dots in the presence of a 
hole or an exciton were calculated based on a microscopic 
theory*2£ Other theoretical studies of quantum dots with 
Mn atoms2i~— were focused on electronic, optical, mag- 
netic or transport properties without any discussion of 
spin lifetimes. Numerous studies of spin relaxation in 
quantum dots were restricted to quantum dots without 
a Mn atorn^ 2 . - — or on diluted magnetic semiconductor 
quantum dots with many Mn atoms 4^ 

In this work, we calculate the relaxation times of Mn 
and electron spin caused by interaction with phonons in 
a singly negatively charged CdTe quantum dot contain- 
ing one Mn atom. Our calculation is based on a micro- 
scopic theory that links the system geometry with spin 
relaxation times. In Sec. |TT] we introduce the theoretical 
model used to describe the system at hand and relevant 
spin relaxation times. In Sec. IIIII we present the results 
and analyze the effects of different electron-phonon in- 
teraction mechanisms, spin-orbit (SO) interaction, mag- 
netic field, temperature, quantum dot dimensions and 
Mn atom position. In Sec. IIVI we compare our results to 
theoretical results for spin relaxation in somewhat simi- 
lar systems and experiments, and analyze the strength of 
other possible spin relaxation mechanisms not included 
in our model. 



II. THEORETICAL MODEL 

In this section, we describe a theoretical model used 
to describe a CdTe quantum dot that contains one extra 
electron and one Mn atom placed in its interior. Such 
a scenario is experimentally realized in a layer of self- 
assembled quantum dots grown at relatively low dot den- 
sity such that interdot interactions are negligible. 
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The Hamiltonian of the system reads 
H = H e i + H m _ e [ + H B , 



(1) 



where H e i is the electronic Hamiltonian, H m _ e i describes 
the interaction of electron with Mn spin and Hb is the 
Zeeman term that describes the interaction of electron 
and Mn spin with external magnetic field. This Hamilto- 
nian acts in the Hilbert space of the system which is given 
as the direct product of electron orbital space, electron 
spin space and Mn spin space. 
The first, electronic term reads 



H e i — ^ Ei. 



at 



(2) 



with c] a and Ci i<T representing electron creation and anni- 
hilation operators. Typical self-assembled quantum dots 
are much larger in the lateral plane (xy plane) than in the 
growth (z) direction, while they can take various shapes 
- lenses, pyramids, truncated pyramids, etc. For this rea- 
son, we adopt the simplest possible quantum dot model 
that captures all essential features of the single-particle 
electronic spectrum - a rectangular box with in-plane di- 
mensions (L x and L y ) much larger than the dimension in 
the z direction (L z ) with infinite potential barriers out- 
side of the region < x < L x , < y < L y , < z < L z . 
The electronic states in the conduction band of semi- 
conductor nanostructures within first several hundreds 
of meV are well described with envelope function effec- 
tive mass Hamiltonian. For our quantum dot model, the 
conduction band electron envelope function is then given 
as 



ip(x, y, z) = N c sin x I sin y ) sin -p- 



n z ir 
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L 

(3) 

where N c = y/8/ (L x L y L z ) is the normalization constant 
and n x , n y and n z are positive integers that represent or- 
bital quantum numbers in each direction. Single particle 
energy Sj ;Cr of an electronic state i with orbital quantum 



numbers 



and 



n z is given as 




(4) 



where m* is the effective mass of the conduction band 
electron in CdTe. 

The second term H m . e i describes the exchange interac- 
tion between an electron and a Mn atom, and it is given 
by the spin impurity model Hamiltonian^ a model well 
established in previous theoretical studies of CdTe quan- 
tum dots with a few Mn atoms2ii22r^5j22i21 



1 



+Z}Ji,tM + + 4 it d u M- . (5) 



In Eq. (refeq:jijr), R is the position of the Mn atom and 
Jy (R) is the electron-Mn spin coupling strength, equal 

to J G V>*(R)V>j(R-)' M z , M+ and M~ are the Mn spin 
operators whose properties are governed by the spin 5/2 
algebra. 

The last term of H a , also known as the Zeeman term, 
describes the interaction of the whole system with an 
external magnetic field B parallel to the z direction. It 
is given as 



H B = -^B9eBS z - \iB9unHM z 



(6) 



where [Ib is the Bohr magneton, and g e (gMn) is the 
electron spin (Mn spin) g factor. S z is the operator of 
the z component of the electron spin. 

In Sec. MI CI we also consider electronic SO cou- 
pling that arises in materials lacking inversion symmetry 
(Dresselhaus SO coupling^) 



H so = 7h • cr, 



(7) 



where as are Pauli matrices, h is the Dresselhaus effec- 
tive magnetic field 

h = [k x (ky — k z ), k y (k z — k x ), k z (k x — k y )], (8) 

with k = — i V and 7 is the Dresselhaus SO coupling 
strength. 

After numerically solving the Hamiltonian eigenvalue 
problem, we obtain the eigenenergies E a and eigenstates 

l*a) 



H \* a ) 
l*a> 



E a \*a), 



(9) 



We use i and S z to denote the orbital and spin state of 
the electron and M z for the Mn spin state. 

Electron-phonon interaction is considered to be the 
main mechanism responsible for the transitions between 
the eigenstates |^ a ) and consequently spin relaxation. 
We will show in Sec. IIII Al that relevant transition ener- 
gies are of the order of meV, which are typical energies of 
acoustic phonons. Because of the high energy of optical 
phonons compared to the acoustic ones, we consider only 
the interaction with acoustic phonons. The Hamiltonian 
of electron-phonon interaction is given a o 7 ' 44 



He-ph — 



E 



(10) 



where w and b are phonon creation and annihilation op- 
erators and q is the phonon wave vector. For acoustic 
phonons, a linear dispersion relation connects the phonon 
wave vector and its energy u) — qv, where v is the sound 
velocity for a particular acoustic phonon branch in a 
given material. 

The scattering matrix depends on the type of 

electron-phonon interaction. For the interaction through 
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deformation potential it is given as^ 
,2 ^ 2 |q| 



\M, 



q,A| 



2Vpv LA ' 



(11) 



and for the interaction through piezoelectric field it can 
be represented as^i 

|Mq ' Al -^n^ (12) 

for longitudinal phonons and 



\M, 



q,A| 



2 2 



a2«5 - 



2 2 

OS 



(3q x q y q z 



(13) 

for transversal ones. For transversal acoustic phonons 
(TA) there are two branches (A=2, 3), whereas for longi- 
tudinal acoustic (LA) phonons there is only one branch 
(A=l). In the preceding expressions D is the acoustic de- 
formation potential, p is the CdTe material mass density, 
V is the volume of the system and £ is given as 



327r 2 e 2 ^ 4 
nVp 



where n is the static dielectric constant and hn is the 
piezoelectric constant. 

For the work presented here, we have used the follow- 
ing parameters: J c — 15eVA 3 , m* — 0.106 m (Ref. I2H) . 

So = -1.67, fl M n = 2.02 (Ref. 25), 7 =11.74eVA 3 
(Ref. Eg), i; L a = 3083m/s, w TA = 1847m/s, 

h u = 3.94 x 10 8 V/m (Ref. EH), D = 5.1 eV, 
p = 4.85 x 10 3 kg/m 3 , k = 9.6 (Ref. E3)- Unless 
otherwise stated, quantum dot dimensions were taken as 
L x = 150 A, L y = 140 A, L z = 30 A. In our calculations, 
we consider first N — 10 electron orbitals, which gives 
the Hilbert space of A^ max = 120 basis vectors. This 
ensures that the calculated quantities have converged to 
their real values. 

Because of the fact that the electron-acoustic phonon 
interaction term is much smaller compared to the rest of 
the Hamiltonian, this term can be treated as a perturba- 
tion responsible for transitions between eigenstates \ty a ) 
of H a . If the system starts in an initial state it will 
make a transition to a final state /) due to electron- 
phonon scattering. The scattering rate for this process is 
determined by Fermi's golden rule 

2tt 



q,A 

' 1 ™q,A + l ± l) 6 ( E f ~ E < ± . (14) 

where Ei and Ef are the energies of the unperturbed 
system in the initial and the final state, is the 

mean number of phonons at a given temperature and 
|( - 4 , /|e iqr |'$i)| is the form factor for electron-phonon in- 
teraction. The (+) sign in Eq. (TTJJ corresponds to the 



process of phonon emission, while the (— ) sign corre- 
sponds to phonon absorption. 

We will show in Sec. IIIII that most of the eigenstates 
of H have a well defined Mn and electron spin because 
one of the cis^m^ coefficients in Eq. © is typically sig- 
nificantly larger than the others. Nevertheless, the re- 
maining coefficients, that correspond to basis states with 
other values of S z and M z , are nonzero. For this rea- 
son, the transitions between the states with different 
values of electron (or Mn) spin are allowed despite the 
fact that the electron-phonon interaction Hamiltonian is 
spin-independent. 

To obtain the average electron spin relaxation time 
one has to consider all possible transitions in the system 
where a particular change of spin occurs. For example, 
the mean relaxation time for electron spin change from 
the initial spin S z — 1/2 to the final spin S z = —1/2 is 
given as 
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Si fi 



(15) 



where the sum over i includes all possible initial states 
with S z = 1/2 and the sum over / includes all possible 
final states with S z = —1/2, with fi being the thermal 
weighting factor of state i at a temperature T, given as 
f i = e X p[-E i /{k B T)]. 



III. RESULTS 
A. Energy spectrum 

The energy spectrum of CdTe quantum dots with a 
single Mn atom has been studied in the pas t 21 ' 30 i 48 and 
here we only review the main features, with an emphasis 
on those that are relevant for our work. 

Relevant energies for our problem are the single- 
particle electron orbital energies Ei, a [Eq. Q], the 
electron-Mn spin exchange interaction energy Jy(R) 
[Eq. ([5])] and the Zeeman splitting energy p, B B [Eq. ©]. 
The separation between the first two orbital energies is 
of the order of 50 meV and is much larger than the ex- 
change interaction and Zeeman splitting energy which are 
of the order of meV. As a consequence, the first twelve 
eigenstates of our system all originate from the ground 
orbital state and are well separated from higher excited 
states, whose average populations are much smaller even 
at room temperature. The dependence of their energies 
on magnetic field is shown in Fig. QJa). 

Looking at the structure of energy levels for electron- 
Mn atom system in Fig. [lja), several distinctive features 
can be noticed. Spin of an electron is 1/2 and that of a 
Mn atom is 5/2. Combined, they will give two possible 
total spin numbers, F = 2 and F = 3, respectively with 
five and seven spin projections along the direction of the 
external magnetic field. When there is no magnetic field, 
total spin is a good quantum number and the presence 
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Figure 1. Magnetic field dependence of the: (a) energies of 
the first twelve energy levels produced by Zeeman splitting 
and electron-Mn exchange interaction; (b) expected value of 
the Mn spin along the direction of an applied magnetic field. 
The system is considered without SO interaction and the Mn 
atom is placed at R = (7.4, 6.9, 1.5) nm. 

of electron-Mn exchange interaction leads to a splitting 
of the twelve-fold degenerate ground level into two new, 
seven-fold and five-fold degenerate ones. In the case of a 
finite magnetic field B, only F z remains a good quantum 
number. The Zeeman term in the Hamiltonian eliminates 
all degeneracies and twelve separate nondegenerate levels 
emerge. 

A typical state of this system is a superposition of all 
basis vectors with the same total spin projection num- 
ber F z . Because electron-Mn spin coupling is relatively 
weak, there are only few dominant states in this linear 
combination. The eigenstate wave function |^) can be 
represented as 

|*) =a\0,+,M-)+P\0 ) -,M+)+..., (16) 

where the basis vectors are ordered by the strength of 
their contribution to the eigenstate. In Eq. (|16[) . |0) de- 
notes the ground orbital state, |+) and — ) are the elec- 
tron spin up and down states, while the Mn spin quantum 
number is related to F z via M±—F z ±1/2, where F z — 
—2, . . . , 2. The next term in Eq. (fTo| depends mostly on 
the Mn-electron coupling strength Joj(R), i.e. on the 
position of the Mn atom. 

For the study of Mn (and electron) spin relaxation, 
the next property that we should turn our attention to is 
the expected value of the Mn spin projection, M z . Its 
dependence on magnetic field for first ten levels with 
F z = —2, . . . , 2 is presented in Fig. W[b). For most val- 
ues of the magnetic field the expected values of M z are 
very close to the corresponding half integer values from 
the interval —5/2 to 5/2, suggesting that Mn spin is well 
defined for a given eigenstate. This corresponds to the 
case where one of the a and (3 coefficients in Eq. (TTB1) is 
much larger than the other. In such a case, the electron 
spin projection S z = F z — M z is also well defined. M z 
and S z cease to be well defined only for certain fields 
and for some states where energy level crossings occur 



(see Fig. [T]). Spin relaxation times calculated within our 
approach should be taken with caution in such cases. 



B. Spin relaxation time 

Since electron-phonon interaction Hamiltonian is in- 
dependent of electron and Mn spin, phonons can only in- 
duce transitions between energy levels with the same to- 
tal spin projection number. Under these conditions, there 
are only five downhill (as well as five uphill) transitions 
allowed. Each of these downhill transitions corresponds 
to Mn spin-flip from M z = F z + 1/2 to M z = F z - 1/2 and 
in the same time to electron spin-flop from S z = —1/2 
to S z = 1/2. 

The transition time between the two states with the 
same F z is therefore also the relaxation time for Mn spin- 
flip from M z = F z + 1/2 to M z = F z - 1/2. On the other 
hand, to obtain the electron spin relaxation time, one has 
to take the average over all possible transitions that lead 
to an electron spin-flip or spin-flop. 

The dependence of Mn and electron spin relaxation 
times on magnetic field in the case of a Mn atom placed 
near the center of the dot [at R = (7.4,6.9, 1.5) nm] at 
room temperature is shown in Figs. [Ha) and[5Jb). 

External field affects the degree of Zeeman splitting, 
which as a consequence determines the energy of the 
phonon through which the system can relax. Since a lin- 
ear dispersion relation connects the phonon energy and 
its wave vector, the external field impact on spin re- 
laxation times comes mostly from scattering matrix ele- 
ments [defined in Eqs. (fTTjl - CTTJl) ]. For relaxation through 
deformation potential, the scattering matrix element is 
~ q, while for piezo-field it is ~ 1/q. Along with the q 2 
factor that comes from the integration over q in Eq. (|14l) , 
this gives an overall ~ 1/q 3 dependence for relaxation 
time through deformation potential and ~ 1/q depen- 
dence for the relaxation time through piezo-field. As a 
consequence, spin relaxation becomes more probable as 
the external field increases, as can be seen from Fig.^a). 
The above mentioned ~ 1/q 3 and ~ 1/q dependences 
are only approximately followed because the form factor 
[Eq. HU] also depends on q but this dependence is rela- 
tively weak in the range of magnetic fields of our interest. 

Fig. [Ha) shows that the results for all possible tran- 
sitions are very similar, ranging from microseconds to 
nanoseconds when the magnetic field varies. This comes 
from the fact that phonon energies for each transition 
are very similar. At zero magnetic field, total spin F (in 
addition to F z ) becomes a good quantum number. Due 
to independence of electron-phonon interaction Hamil- 
tonian on spin, the transitions between the states with 
cither different F or F z become forbidden. This leads to 
infinite relaxation times at B = in our model. 

Besides the Mn spin, the electron spin relaxation also 
occurs. As mentioned above, a consequence of the spin 
conserving Hamiltonian is the connection between these 
two. Each time a Mn atom makes a flip to a neighboring 
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Figure 2. Spin relaxation times for a Mn atom placed 
at R = (7.4,6.9,1.5) nm at room temperature (T=295K), 
caused by deformation potential interaction with LA phonons 
(dl), piezoelectric potential interaction with LA (pi) and TA 
phonons (pt): (a) Five transitions corresponding to Mn atom 
spin-flip, (b) Electron spin relaxation time obtained by av- 
eraging. Empty symbols represent the data for spin-down to 
spin-up process, while the filled ones represent the opposite, 
(c) The same as (b) but for the calculation performed using 
an inaccurate approximation. 



spin state, the electron makes a Hop and therefore the 
electron spin relaxation time is obtained by averaging 
over all possible transitions. The similarity between these 
two can be seen from Fig. HJb). There are two types of 
processes for electrons: from spin-down to spin-up state, 
and the opposite one. In general there is a difference 
between the relaxation times for these two but only in 
the case of sufficiently low temperatures. 

As mentioned earlier, one, or at most two, coefficients 
in Eq. (|16[) give the dominant contribution to an eigen- 
state. Both of these coefficients correspond to the ground 
orbital state. It is therefore very tempting to introduce 
an approximation in which we would reduce the Hilbert 
space of the system to basis states that originate from 
orbital ground state only. However, this is not appro- 
priate since such an approximation would yield infinite 
relaxation times. A less drastic approximation, where the 
eigenstates are calculated accurately, but only first two 
terms in Eq. (fTOl) are kept for the calculation of transition 
rates, is also highly inaccurate. As shown in Fig. [2fc), 
it gives unrealistically large spin relaxation times. The 
results obtained from these tempting, but inaccurate ap- 
proximations, demonstrate the necessity of including a 



larger number of orbital states in the calculation and the 
need for a numerical, rather than analytical, approach to 
the problem. 



C. The role of SO interaction 

The electron (Mn) spin relaxation processes caused by 
spin-independent electron-phonon interaction are only 
possible due to the presence of terms in H Q that mix 
the states of different spin. Such a term in H a is the 
electron-Mn exchange interaction H m - e U which mixes 
both the states of different electron spin and the states 
of different Mn spin. 

Another term which leads to mixing of the states of 
opposite electron spin is the SO interaction [Eq. ([?])]■ 
The results presented in Sec. IIIIBI were obtained while 
neglecting this term and the goal here is to assess the 
accuracy of such an approximation. 

In the presence of SO interaction, F z is no longer a 
good quantum number. However, it turns out that the 
mean value of F z for an eigenstate is very close to the 
value of F z in the absence of SO interaction. Therefore, 
the eigenstates of the Hamiltonian can still be labeled by 
the mean value of F z for a state. 

In Figs. G^a) and (c) we compare the electron spin re- 
laxation times in the presence and absence of SO interac- 
tion. The values appear to be rather similar which sug- 
gests a weak effect of SO interaction on spin relaxation 
times. 

A qualitative difference that SO interaction introduces 
is that it allows for transitions between the states with 
different values of F z , since F z is no longer a good quan- 
tum number. The calculated transition rates for such 
transitions from the lowest state with F z = 2 to other 
states are shown in Fig.[3Jb). These transitions appear to 
be much weaker than the transitions between the states 
with the same F z . 

As a conclusion to this section, we may say that the 
changes in spin relaxation times due to SO interaction 
are very small and for all practical purposes SO interac- 
tion can be neglected. Therefore, the rest of the results 
presented in this paper will not include SO interaction. 



D. Mn position and temperature 

The position of Mn atom determines the electron-Mn 
exchange coupling constants Jy (R) [Eq. which is 
the only place where it appears in the Hamiltonian. Since 
this is the only term in the Hamiltonian that mixes the 
states of different electron or Mn spin, one may expect 
that it has a significant effect on spin relaxation times. 
To understand the role of the Mn atom position on spin 
relaxation times, we have performed a calculation for two 
positions of Mn atom. Position 1 is near the center of the 
dot at R = (7.4, 6.9, 1.5) nm, while position 2 was chosen 
to maximize the Joi(R-) coupling constant and is given 
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Figure 3. The influence of SO interaction on spin relaxation times for a Mn atom placed at R = (7.4, 6.9, 1.5) nm. Magnetic 
field dependence of relaxation time caused by deformation potential interaction with LA phonons for: (a) the transition between 
two levels with the same mean value of F z in the presence (full symbols) and absence (empty symbols) of SO interaction; (b) 
the transitions between the level (F = 3, F z = 2) and other levels. Full line represents the data for the transition to (F = 2, 
F z =2), while dashed, dotted and dash-dotted lines correspond to transitions to (F = 3, F z = 1), (F — 2, F z — 1) and (F — 3, 
F z = 3) respectively. The transitions to other levels that are not shown have relaxation times longer than 1 ms. The levels were 
labeled according to the quantum numbers F and F z which are good quantum numbers at zero magnetic field in the absence of 
SO interaction, (c) Magnetic field dependence of average electron spin-up to spin-down relaxation time in the presence (filled 
symbols) and absence (empty symbols) of SO interaction. The labels dl, pi and pt have the same meaning as in Fig. [2] 



as R = (4.5,7,1.5) nm. The obtained spin relaxation 
times are shown in Fig. |4j At position 1, the Joi, J02 
and J03 coupling constants are nearly zero. As a conse- 
quence, the states of interest have the first largest contri- 
bution in the linear combination [Eq. (Tl6l) ] other than a 
and /3 from coefficients corresponding to orbital state |4) 
\ c 1s m coefficients in the notation of Eq. ©]. On the 
other hand, at position 2, the first largest contribution 
other than a and j3 comes from c1 S M . Since state |1) 
has a lower energy than state |4), the coefficients c^ s M 
for position 2 are larger than c|c M coefficients in the 
case of position 1. As a consequence, the relevant states 
at position 2 exhibit a stronger mixing of basis states of 
different spin, which translates into faster spin relaxation 
at position 2, as can be seen from Fig. @] 
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Figure 4. The influence of Mn atom position and temperature 
on electron spin relaxation times. The dependence of spin 
relaxation time for the transition from S z = 1/2 to S z = —1/2 
on magnetic field for different temperatures T and Mn ion 
positions R - case 1: T = 295 K, R = (7.4,6.9,1.5) nm; 
case 2: T = 3 K, R = (7.4, 6.9, 1.5) nm; case 3: T = 295 K, 
R= (4.5,7.0,1.5) nm. 



Spin relaxation times at temperatures of 3 K and 295 K 
are also shown in Fig. |4j When an individual transi- 
tion is concerned, the temperature appears in our theory 
only through the phonon occupation number and con- 
sequently its effect is easily predictable - higher tem- 
peratures lead to shorter spin relaxation times. When 
average spin relaxation times are concerned, the tem- 
perature appears in the theory also through the thermal 
weighting factors [Eq. (fT5"j) ]. This becomes especially im- 
portant at higher magnetic fields, where the state with 
F z — 3 becomes the ground state. The transition prob- 
ability from this state to other states is small but this 
state has a high weighting factor. The average transition 
rate is then significantly different than the transition rate 
for individual transitions between states with the same 
F z (F z = —2, ... ,2). As a consequence of all the men- 
tioned temperature effects, cooling down the system from 
room temperature to cryogenic temperature leads to an 
increase in spin relaxation times by three orders of mag- 
nitude. 



E. Quantum dot dimensions 

The dependence of spin relaxation times on quantum 
dot dimension L z is shown in Fig. [5ja). The relative 
position of an Mn atom inside the dot is kept during 
this change of dimensions. The electron orbital states 
[Eq. (|3J)] in the range of energies that is of interest here 
all have quantum number n z = 1. As a consequence, 
their energies do not change when L z is changed (upto 
an irrelevant constant). Therefore, the change in L z af- 
fects only the electron wave functions and consequently 
the electron-Mn exchange coupling constants Jy(R). 
When the dot dimensions and the confinement volume 
increases, the probability of finding an electron near the 
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Mn atom, due to normalization of the wave function, 
decreases. As a consequence, Jjj (R) constants decrease. 
This leads to weaker mixing of states of different spin and 
therefore spin relaxation times increase, when quantum 
dot dimensions increase, as can be seen in Fig. Ufa). 




L Z (A) L x (A) 

Figure 5. The influence of quantum dot dimensions on elec- 
tron spin relaxation times at a temperature of T = 295 K 
and magnetic field of B = 5 T. Filled symbols represent 
the data for spin-up to spin-down process, while the empty 
ones represent the opposite. The labels dl, pi and pt have 
the same meaning as in Fig. [2] The Mn atom is placed at 
R = (4.5,7.0, 1.5) nm for the default dot dimensions and its 
relative position is kept constant when dot dimensions are 
varied, (a) The dependence of electron spin relaxation time 
on the dot dimension in z direction; (b) The dependence of 
electron spin relaxation time on the dot size in xy plane - the 
Lx/Ly ratio and L z are kept constant, while L x is varied. 

For the same reasons [weaker Jij(R) for larger dots], 
one may expect that the relaxation times will increase 
when the dot dimensions in xy plane increase. How- 
ever, the trend obtained in Fig. EJb) is somewhat differ- 
ent. The reason is that when L x and L y increase, the 
distance between electron orbital energies [Eq. (gj)] de- 
creases. This leads to stronger contributions from basis 
states that originate from electron orbitals states other 
than |0) [the missing terms in Eq. (|16[) ]. i.e. to stronger 
mixing of states of different spin in the expansion from 
Eq. (|16|) . which leads to shorter spin relaxation times. 
To summarize, an increase of dot dimensions in xy plane 
leads on the one hand to a decrease of Jy (R) and on the 
other hand to an increase in mixing in Eq. (|16|) . As a 
consequence of these two opposite trends, one obtains a 
nonmonotonous dependence of spin relaxation times on 



IV. DISCUSSION 

In this section, we discuss the relevance of other possi- 
ble spin relaxation mechanisms which were not included 
in our model and compare our results to other relevant 
results from the literature. 

In this work, we have only considered the Dresselhaus 



SO coupling which is a consequence of bulk inversion 
asymmetry. Realistic quantum dots also exhibit Rashba 
SO coupling*- as a consequence of structural inversion 
asymmetry. Our model quantum dot is symmetric, there- 
fore Rashba SO coupling is not present. However, real- 
istic quantum dots certainly exhibit a certain degree of 
asymmetry and consequently the Rashba SO coupling. 
The interplay of Dresselhaus and Rashba SO coupling 
in quantum dots has been studied in Ref. [H. It was 
found that both of these are of similar strength and cause 
similar spin relaxation times. Moreover, in Ref. |32| the 
authors found that Dresselhaus SO coupling provides a 
bigger admixture of states of different spin. Since we have 
found that Dresselhaus SO coupling has practically no ef- 
fect on spin relaxation rates in our system, it is expected 
that the same will be the case for Rashba SO coupling. 

Spin relaxation by direct electron spin-phonon cou- 
pling in GaAs quantum dots was also considered in 
Refs. [32| and HO- A conclusion was reached that it is by 
far less effective than the relaxation caused by SO inter- 
action induced mixing and electron-phonon interaction. 

Spin relaxation due to hyperfine interaction between 
the electron s pin and the spin of the nuclei was investi- 
gated in Ref. l33l The hyperfine interaction was found 
to be the dominant mechanism responsible for spin re- 
laxation rates at magnetic fields which are low enough. 
Since our model gives relaxation rates that tend to zero at 
vanishing magnetic field, one may expect that hyperfine 
interaction will become relevant at low magnetic fields in 
our system too. 

Direct coupling of Mn spin to phonons was considered 
in Ref. [2^ as a potential mechanism of Mn spin relax- 
ation. The corresponding relaxation times were found to 
be at least of the order of 1 ms, being much longer than 
the ones originating from other mechanisms in our work. 

Based on the discussion above, we may say that the 
spin relaxation mechanism considered in our work is cer- 
tainly the most relevant mechanism in the system stud- 
ied. 

Next, we compare the trends that we obtain for the 
dependence of relaxation times on various parameters to 
the ones obtained in quantum dots without Mn atoms. 
In Ref. [36|, spin relaxation in GaAs quantum dots caused 
by the presence of SO interaction and electron-acoustic 
phonon interaction was investigated. In that situation, 
spin relaxation takes place due to the transition between 
the first two Zeeman sublevels of opposite spin. The only 
spin mixing mechanism in such a system is the SO in- 
teraction. Spin relaxation times were found to decrease 
with the increase in magnetic field and were found to 
be mainly determined by the value of the phonon wave 
vector responsible for the transition. These conclusions, 
as one might have expected, are the same in our work. 
Interestingly, it was found on the other hand in Ref. [36| 
that the dominant electron-phonon interaction mecha- 
nism that causes spin relaxation is piezoelectric interac- 
tion with TA phonons. In our system, piezoelectric in- 
teraction with TA phonons is comparable to deformation 
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potential interaction with LA phonons at low magnetic 
fields, whereas at higher magnetic fields deformation po- 
tential interaction with LA phonons becomes the dom- 
inant relaxation mechanism. When the dependence of 
relaxation times on dot dimensions in xy plane is con- 
cerned, it was found in Ref. that these decrease when 
dot dimensions increase. The presence of exchange inter- 
action between the electron and Mn spin, as the domi- 
nant spin mixing mechanism in our system, causes a dif- 
ferent trend in our case, as seen in Sec. lIIIEl It was found 
that spin relaxation times increase with the increase of 
dot dimensions in the z direction both in our work and 
in Ref. 36, albeit for a different reason. In our case, this 
is caused by the decrease of electron-Mn spin exchange 
coupling constant when the dot volume increases (see 
Sec. IIII D|) . while in the work of Ref. |36j this is caused by 
the decrease of SO interaction. 



Finally, we compare the electron spin relaxation times 
that we obtained to the limits on their values in- 
ferred from Mn spin optical orientation experiments. 17 
In Ref. |H, a theoretical analysis of Mn spin optical ori- 
entation experiments^ was performed to understand the 
physical origin of the observed Mn spin orientation. It 
was found that optical orientation that occurs on the 
~10ns timescalei 7 - can be explained if hole spin relax- 
ation times of the order of ~10ns are assumed. Similar 
hole spin relaxation times were obtained from a micro- 
scopic theory in Ref. H^J. In the same time, the electron 
spin relaxation time needs to be longer than that. 19 Our 
theoretical results at a low temperature and low magnetic 
field (B < 1 T) (Fig. 3|) indicate electron spin relaxation 
times longer than tens of microseconds, which is fully 
consistent with the conclusion obtained in Ref. [l9|. 



V. CONCLUSION 

In conclusion, we found that the interactions responsi- 
ble for electron and Mn spin relaxation in quantum dots 
doped with a single Mn atom are the electron-Mn spin 
exchange interaction and the electron-phonon interac- 
tion. The former provides mixing of the states with dif- 
ferent spin within an eigenstate and allows for the spin- 
flip transition caused by the latter. SO interaction has 
a negligible effect on spin relaxation times, in contrast 
to conventional quantum dots where SO interaction is 
the only mechanism which allows for mixing of states 
of different spin within an eigenstate and the relaxation 
through the interaction with phonons. Among the differ- 
ent electron-phonon interaction mechanisms, electron in- 
teraction with LA phonons turns out to be the dominant 
mechanism responsible for spin relaxation. Spin relax- 
ation times decrease with the increase in magnetic field. 
This dependence is mainly determined by an increase in 
energy level splitting and consequently the increase in the 
wave vector of the phonon responsible for the transition. 
The position of Mn ion within the dot determines the 
strength of electron-Mn spin exchange interaction and 
therefore significantly alters the spin relaxation times. 
We find that spin relaxation times in our electron-Mn 
spin system are longer than in similar hole-Mn spin or 
exciton-Mn spin systems. This suggests that potential 
spintronic or quantum computing devices based on the 
interaction between charge carriers and Mn spin should 
use electrons as charge carriers. 
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